(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 

International Bureau 




llllllllllllllllllllllllllllllllllllllllllllllli 



(43) International Publication Date (10) International Publication Number 

16 May 2002 (16.05.2002) PCT WO 02/39136 A2 



(51) International Patent Classification 7 : G01S 

(21) International Application Number: PCT/US01/32163 

(22) International Filing Date: 16 October 2001 (16.10.2001) 

(25) Filing Language: English 

(26) Publication Language: English 

(30) Priority Data: 

09/691,690 18 October 2000 (18.10.2000) US 

(71) Applicant: HONEYWELL INTERNATIONAL INC. 

[US/US]; 101 Columbia Avenue, P.O. Box 2245, Morris- 
town, NJ 07960 (US). 

(72) Inventor: BRENNER, Mats, A.; 4275 Deerwood Lane, 
Plymouth, MN 55441 (US). 



(74) Agents: CRISS, Roger, H. et ah; Honeywell International 
Inc., 101 Columbia Avenue, P.O. Box 2245, Morristown, 
NJ 07960 (US). 

(81) Designated States (national): CA, JP. 

(84) Designated States (regional): European patent (AT, BE, 
CH, CY, DE, DK, ES, H, FR, GB, GR, IE, IT, LU, MC, 
NT,, PT, SE, TR). 

Published: 

without international search report and to be republished 
upon receipt of that report 

For two-letter codes and other abbreviations, refer to the ''Guid- 
ance Notes on Codes and Abbreviations " appearing at the begin- 
ning of each regular issue, of the PCT Gazette. 



(54) Title: APPARATUS FOR NAVIGATION SATELLITE SIGNAL QUALITY MONITORING 



IP 



•El 



IL1 



< 



IE 2 



IL 2 



(57) Abstract: An apparatus for the detection of posi- 
tioning system satellite signal distortions includes a cor- 
relator that determines a plurality of correlation measure- 
ments at points along a correlation curve. The correlation 
measurements are based upon a correlation between a re- 
ceived satellite signal and a reference. A signal distortion 
detector determines differences between the correlation 
measurements along the correlation curve and detects a 
signal distortion from the differences. 
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APPARATUS FOR NAVIGATION SATELLITE 
SIGNAL QUALITY MONITORING 

5 

Technical Field of the Invention 

The present invention relates generally to satellite based 
positioning systems such as the Global Positioning System (GPS) and more 
particularly to the monitoring of the quality of the signals transmitted by 
10 satellites in a satellite based positioning system. 

Background of the Invention and Prior Art 

A satellite based positioning system is used to determine a 
position of a receiver and typically includes satellite control facilities, a plurality 

15 of satellites, the receiver, and one or more local or regional ground stations. 
Each of the satellites transmits a signal that contains a code and certain 
prescribed information useful to the receiver in determining its position. The 
receiver synchronizes itself to the codes of at least four satellites and uses the 
information in the signals from these satellites in order to perform a 

20 triangulation like procedure so as to determine its coordinates and time offset 
with respect to a reference, such as the center of the Earth and the GPS 
standard time. 

The receiver is not constrained to a specific location and, 
therefore, represents a variable position. Indeed, the purpose of the satellite 

25 based positioning system is to make it possible for the receiver to determine its 
position regardless of the iocation of the receiver. On the other hand, the local 
or regional ground station is in a fixed location and is used to monitor the 
signals transmitted by the satellites. The signals transmitted by the satellites 
can be adversely affected, for example, by atmospheric conditions which can 

30 lead to improper position determinations by the receiver. The ground station, 
therefore, notifies the receiver of any necessary signal corrections to allow the 
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receiver to make more accurate position calculations. This arrangement is 
referred to as differential positioning. 

The ground station of the present invention also monitors the 
signals transmitted by the satellites in order to detect faults within the satellites. 
5 For GPS, these faults are specified by the FAA who imposes stringent 

requirements to protect users against positioning system signal faults. A set of 
test waveforms has been chosen by the FAA to represent at least some of the 
more egregious faults. These waveforms are used for certification testing of 
the ground station equipment. 

10 The prior art determines faults by comparing conventional code 

tracking discriminators at different correlator spacings. As shown in Figure 1 , a 
correlation curve is established by correlating the code received from a satellite 
with a suite of code references which are time shifted replicas of the code 
transmitted by that satellite. For example, seven correlation measurements 

15 may be calculated as shown in Figure 1 . The in-phase measurement IP 
represents the amount of correlation between the received code and a 
reference code that has a zero time shift with respect to the received code (this 
measurement is referred to as punctual). The in-phase measurement IE 1 
represents the amount of correlation between the received code and a 

20 reference code that has a first predetermined time shift so that it is early with 
respect to the received code. The in-phase measurement ILi represents the 
amount of correlation between the received code and a reference code that has 
a second predetermined time shift so that it is late with respect to the received 
code. Similarly, the in-phase measurement IE 2 is derived using a third 

25 predetermined time shift, the in-phase measurement IL 2 is derived using a 

fourth predetermined time shift, the in-phase measurement IE 3 is derived using 
a fifth predetermined time shift, and the in-phase measurement IL3 is derived 
using a sixth predetermined time shift. The magnitude of the first 
predetermined time shift may be equal to the magnitude of the second 

30 predetermined time shift, the magnitude of the third predetermined time shift 
may be equal to the magnitude of the fourth predetermined time shift, and the 
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magnitude of the fifth predetermined time shift may be equal to the magnitude 
of the sixth predetermined time shift. It is assumed that all measurements are 
normalized such that the measured correlation is a function of the time shifts 
only and not the absolute power of the received satellite signal. 
5 First, second, and third discriminators are then formed 

according to the following equations: 

dsubl ~ = ~(lLsub 1 ~-~IEsub1)IP 

10 d sub 2 ~ = ~~ (IL sub 2 ~ - ~ IE sub 2) IP 

d sub 3 ~ = ~ (IL sub 3 ~ -~ IE sub 3) IP 

These discriminators are thereafter compared to each other through the 
15 formation of quantities di j2 , di t 3, and d 2|3 according to the following equations: 

d sub {1 ,2} ~ = ~ LINE d sub 1 ~ - ~ d sub 2 LINE 
d sub {1 ,3} ~ = ~ LINE d sub 1 ~ - ~ d sub 3 LINE 

20 

d sub {2,3} ~ = ~ LINE d sub 2 ~ - ~ d sub 3 LINE 



The quantities d^ ch^, and d2,3 are compared to corresponding thresholds 
Di,2, D 1j3 , and D 2 ,3 such that, if the first discriminator di >2 exceeds D 1>2 , if the 
25 second discriminator di )3 exceeds Di ^ or if the third discriminator d 2j 3 exceeds 
D 2>3 , a fault is assumed to exist During normal operation of the global 
positioning system, this test is performed on the signals received from each of 
the satellites. During certification, a test is to be performed using each. of the 
test waveforms chosen by the FAA in order to prove that fault detection occurs. 
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At least one of the problems with this method is that it is 
requires six correlators in order to determine the three quantities d^, di,3> and 
d2,3 which is too much hardware for the amount of useful data being provided. 

It is also known for ground stations to determine faults by 
5 scanning the whole correlation peak (i.e., the portion of the correlation curve 
around the punctual in-phase measurement IP) in order to determine whether 
the peak varies from some prescribed norm by a predetermined amount. 
However, this fault detection arrangement requires a substantial amount of 
computing power and it lacks accuracy. 
10 A third method in the prior art uses the following ratios between 

the measurements IE3, IE2, IE1, IL1, IL2, and IL3: 

r sub { E3,E2} ~ = ~ {IE3} over {IE2} 

15 r sub { E3,E1} ~ = ~ {1E3} over {IE1} 

r sub { E3,L1} ~ = ~ {1E3} over {IL1} 



20 



Each of these ratios is compared to a corresponding predetermined value. 

The present invention is directed to an arrangement which 
25 overcomes one or more problems of the prior art. 

Summary of the Invention 

In accordance with one aspect of the present invention, an 
apparatus for the detection of positioning system satellite signal faults 
30 comprises a correlator and a fault detector. The correlator determines a 

plurality of correlation measurements at points along a correlation curve, and 
each correlation measurement is based upon a correlation between a received 
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sateliite signal and a reference. The fault detector determines differences 
between the correlation measurements along the correlation curve and detects 
a fault from the differences. 

In accordance with another aspect of the present invention, a 
5 method of detecting faults affecting a signal transmitted by a positioning system 
satellite comprises: correlating the transmitted signal with a first reference in 
order to determine a first correlation measurement at a first point along a 
correlation curve; correlating the transmitted signal with a second reference in: 
order to determine a second correlation measurement at a second point along 

10 the correlation curve; correlating the transmitted signal with a third 'reference in .i 
order to determine a third correlation measurement at a third point along the 
correlation curve; determining a first difference from the first and second J: ; .:■ 
correlation measurements; determining a second difference from the second 
and third correlation measurements; directly comparing the first difference to a 

15 first threshold; directly comparing the second difference to a second threshold; i 
and, detecting a fault in the satellite based upon the comparisons of . the first * 
and second differences to the first and second thresholds. 

In accordance with still another aspect of the present invention, 
a method of detecting faults affecting a signal transmitted by a positioning 

20 system satellite comprises: correlating the transmitted signal with references in 
order to determine a plurality of correlation measurements at corresponding 
points along a correlation curve; determining a single value from n pairs of the 
correlation measurements, wherein n > 2; comparing the single value to a 
threshold; and, detecting a fault in the satellite based upon the comparison. 

25 

Brief Description of the Drawings 

These and other features and advantages will become more 
apparent from a detailed consideration of the invention when taken in 
conjunction with the drawings in which: 
30 Figure 1 is a waveform showing a correlation diagram useful in 

explaining prior art fault detection as implemented in a ground station in a 
global positioning system; 
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Figure 2 is a schematic diagram of a portion of a ground station 
receiver pertinent to the present invention; and, 

Figure 3 is a waveform showing a correlation diagram useful in 
explaining fault detection as implemented by a ground station in a global 
5 positioning system in accordance with the present invention. 



Detailed Description 

A portion of a '-ground -station 1.0 pertinent to the present 
nvention is shown in Figure 2;cTT1e ground station has correlators 2-E,vi,- i i s.; 

lb 1 H2-E3, 12-E2, 12-Ei, 12-P; "12-^2-1:2! 12^3, . . ./12-Ln, where n'4-Ms 12-?- 
greater than two, and where n is the number of late correlation measurements 
and m is the number of early correlation measurements to be usedcin o -? :c, 
determining a fault. The correlator 12-P correlates the usual code in the • ; 
received signal with a reference 14-P to produce a punctual correlation output 

15 IP, the correlator 12-L-t, correlates the code in the received signal with a a 
reference 14-L-j to produce a late correlation output IL1, the correlator 12-L 2 
correlates the code in the received signal with a reference 14-L ? to produce a 
late correlation output IL2, the correlator I2-L3 correlates the code ih the \ 
received signal with a reference I4-L3 to produce a late correlation output I L3; . 

20 . and the correlator 12-L n correlates the code in the received signal with a 
reference 14-L n to produce a late correlation output IL n . «■ 

In addition, a correlator 12-Ei correlates the code in the 
received signal with a reference 14-Ei to produce an early correlation output 
IEi, a correlator 12-E 2 correlates the code in the received signal with a 

25 reference 14-E 2 to produce an early correlation output IE 2 , a correlator 12-E 3 
correlates the code in the received signal with a reference 14-E 3 to produce an 
early correlation output IE 3 , . . and a correlator 12-E m correlates the code in 
the received signal with a reference 14-E m to produce an early correlation 
output IE m . 
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The ground station 10 has a processor 16 which uses the 
punctual and late correlation outputs IP, ILi, IL 2 , IL3, . . IL n as disclosed 
hereinafter in order to determined whether a fault exists. Alternatively or 
additionally, the processor 16 can use the early correlation outputs IE 1? IE 2 , IE 3 , 

5 IE m as disclosed hereinafter in order to determined whether a fault exists. 

In order to generate the punctual correlation output IP, the 
processor 16 shifts the reference 14-P, which may be a replica of the code 
contained in the received signal, until an optimum correlation is obtained. The 
processor 16 then controls the reference 14-Li so that the reference 14-Li is a 

0 replica of the reference '14-P and so that the reference 14-L-i is time shifted with 
respect to the reference 14-P by a first predetermined amount of time. 
Accordingly, the correlator 12-L-| produces the late correlation output ILj. The 
processor 16 also controls the reference 14-L 2 so that the reference 14-L 2 is a 
replica of the reference 14-P and so that the reference 14-L 2 is time shifted with 

5 respect to the reference 14-P by a second predetermined amount of time, 
where the second predetermined amount of time is greater than the first 
predetermined amount of time. Accordingly, the correlator 12-L 2 produces the 
late correlation output IL 2 . Similarly, the processor 16 controls the reference 
14-L3 so that the reference 14-L3 is a replica of the reference 14-P and so that 

0 the reference M-L3 is time shifted with respect to the reference 14-P by a third 
predetermined amount of time, where the third predetermined amount of time is 
greater than the first and second predetermined amounts of time. Accordingly, 
the correlator 12-L 3 produces the late correlation output IL 3 . The remaining late 
correlation outputs up to IL n are generated in a like manner. The first, second, 

:5 third, etc. predetermined amounts of time are all chosen so that the late 

correlation outputs I L-i through IL n are all on the downward or late slope of the 
correlation curve as shown in Figure 3. 

Additionally or alternatively, the correlators 12-E-j, 12-E 2 , 12-E 3 , 
. . ., 12-Em nnay be positioned so as to generate the early correlation outputs 
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IEi, IE2, IE3, . . IE m . Also, quadrature phase correlation outputs QE m , . . 
QE1, QP, QL«t . ., QL n may be generated by correlating the code in the 
received signal to a time shifted quadrature form of the reference 14-P. In 
accordance with this latter alternative, each measurement used to generate a 
5 fault indication may be formed as an RMS (Root Mean Square) value of the 
corresponding in phase and quadrature phase measurements. 

The set IE m , . . IE 3 , !E 2 , lE^ IP, IL 1f IL 2 , IL 3 , . . lL n may be 
denoted as l m , . . ., L 3 , L 2 , L*i, lo, U, '2, h, ■ ■ - Jn and the following corresponding 
set of RMS values 

10 

SQRT { IE sub m sup 2 ~ + ~ QE sub m sup 2}, ~ ~ ~ , SQRT { IE sub 1 sup 2 
~ + ~ QE sub 1 sup 2}, ~ SQRT { IP sup 2 ~ + ~ QP sup 2}, ~ ~ ~ , -SQRT { IL 

sub n sup 2 ~ + ~ QL sub n sup 2} 

15 may be denoted as R m , . . ., R_ 3 , R_ 2 , R-i, Ro, R1, R2» R3, ■ ■ Rn- 

If early as well as late correlation outputs are to be used for fault 
detection, the processor 16 processes the early correlation outputs IE m through 
IE-i, the punctual correlation output IP, and/or the late correlation outputs IL 1 
through lL n so as to derive one or more measured differences dy. These 

20 measured differences dy are generated in accordance with the following 
equations: 

d sub {i,j} - = - | sub i I sub j 

25 or 

d sub {i,j} ~ = ~ R sub i ~ - ~ R sub j 

where i = -m, . . n and j = -m, . . n, and where the negative sign indicates 
30 measurements on the early slope and the positive sign indicates 
measurements on the late slope of the correlation curve. 
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At this point, it is possible to subtract the expected difference 
from all or a subset of these measured differences dj j and to compare the 
resulting difference deviations to corresponding thresholds in order to 
determine the existence of a fault. For example, assuming that all of these 
5 difference deviations are used, then these difference deviations may be 
compared to corresponding thresholds in accordance with the following 
equation: 

LINE d sub {i j} ~ - ~ Ed sub {ij} LINE - > - D sub {i j} (: 

10 

where Edg is the difference that is expected for each corresponding measured 
difference djj when there is no fault. 

In some cases, the measured differences djj may be affected by 
thermal and multipath noise which could lead to false detection of faults, 

15 depending upon the sensitivity of the fault detection apparatus, i.e., the 

magnitudes of the thresholds Djj. Accordingly, in these cases, a fault could be 
detected when no fault is in fact present, or a fault which is present might not 
be detected at all. 

The thermal noise content in djj can be determined as a 

20 function of the delay hjj between the reference codes 14-E m , . . ., 14-E3, 14-E 2 , 
14-E-i, 14-P, 14-Li, 14~L 2 , 14-L 3 , . . 14-L n . The delay hjj is the delay between 
the two references that are correlated with the received signal to produce Ij and 
Ij. Typically, hjj = 0.025 to 0.05 chip, but may vary from this range. The 
thermal noise th1 in djj depends on the signal to noise ratio and the standard 

25 deviation (1-sigma) of th1 and is given by the following equation: 

sigma sub {th1} ( ij ) ~ = ~ 293SQRT { { h sub {ij} B } over { S/No } } (4) 

where B is the two-sided bandwidth of the noise. In addition, there is another 
30 contribution, th2, to the thermal noise due to the variation of the punctual 
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reference (i.e., the reference 14-P). Accordingly, the total thermal noise is th = 
th1 + th2. The multipath noise mp depends on the antenna gain pattern and its 
overbounding 1-sigma ^mp(U)( v ert 50 {stack {vert -50 {..} # e}} vert -50 ) is 
expressed as a function of satellite elevation vert 50 {stack {vert -50 {..} # e}}. 
5 The statistical properties of th and mp are usually identified at installation of the 
ground station and the statistical information is parameterized and are 
thereafter stored in memory. 

One way to minimize any adverse effects of thermal and 
multipath noise is to make a plurality of measurements for each of the 

10 measured differences djj that are used in the detection of faults. Then, the 
measurements for each of the measured differences djj may be averaged or 
filtered. Because the thermal noise and some of the multipath noise are not 
particularly correlated from one measurement to the next, averaging will tend to 
reduce the effects of thermal and multipath noise. 

15 As an example, let it be assumed that the punctual correlation 

output IP and the late correlation outputs IL-j and IL2 are used to detect faults. 
Accordingly, the following measured differences are determined: d 0 ,i = IP - IL-i; 
do,2 = IP - IL-2; and, di^ = IL1 - IL-2. In order to reduce the effects of thermal 
and multipath noise, however, plural calculations of the measured difference 

20 d 0 ,i are made based upon plural correlation measurements resulting in plural 
punctual correlation outputs IP and plural late correlation outputs IL1. All such 
calculations of the measured difference do,i are then averaged. Similarly, 
plural calculations of the measured difference do,2 are made based upon the 
plural correlation measurements resulting in plural punctual correlation outputs 

25 IP and plural late correlation outputs IL 2 . As before, all such calculations of the 
measured difference do,2 are averaged. Likewise, plural calculations of the 
measured difference d 1>2 are made based upon the plural correlation 
measurements resulting in the plural late correlation outputs IL-j and plural late 
correlation outputs IL 2 . Again, all such calculations of the measured difference 
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di j2 are averaged. These averages may then be compared to their 
corresponding thresholds D 0j i, D 0 ,2> and D 1j2 in order to determine the 
existence of a fault. 

Another way to reduce the effect of thermal and multipath noise 
5 is to suitably filter the measured differences dg, or the punctual correlation 

output IP, the late correlation outputs II_i through IL n , and the early correlation 
outputs IEi through IE m , such as with a low pass filter. 

Still another way to reduce the effect of thermal and multipath 
noise is by implementing the following procedure. In describing this procedure, 
10 it is useful to define a covariance matrix P in accordance with the following 
equation: 

P ~ = ~ E [ ( underline d underline m ) ( underline d ~ - ~ underline m ) sup 

T] 

15 

where the underlines indicate vectors, where E[A] is the statistical expectation 
of A, where the vector m is the mean value of the vector d, and where the 
vector d is determined in accordance with the following equation: 

20 {underline d} sup T ~ = ~ ( d sub 1 , d sub 2, d sub 3, d sub 4, . . d sub N) (6) 

where dk = Ik - 'k-1 - Ed k for k = -m, . . ., n-1 or where d k = Rk - Rk-1 - Ed k for k = 
-m, . . n-1 assuming N + 1 correlation measurements such as l m , . . L3, L 2 , I. 
1, lo> h, l2> '3, ■ ■ •> 'n- An upper triangular matrix U and a diagonal matrix D are 
25 determined according to the following equation: 

P~ = ~UDUsupT 

where P is the covariance matrix given by equation (6). With the covariance 
30 matrix P known from equation (6), the upper triangular matrix U and the 
diagonal matrix D can be determined, for example, by using Cholesky 
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factorization. Thus, the following relationship may be defined in accordance 
with the following equation: 

stack { vert -50 {\~} # d # vert 50 {-} } ~ = ~ U .sup {-1} ( underline d - - ~ 
5 underline m) 

where stack { vert -50 {\~} # d # vert 50 {-} } is a vector representing the 
decorrelated deviations generating the vector <± Equation (9) can be re-written 
according to the following equation: 

10 

underline d ~ = ~ U ~ stack { vert -50 {\~} # d # vert 50 {-}}- + - underline m(9) 

Then, combining equations (6) and (10) produces the following equation: 

15 P~ = ~E~[U~ stack { vert -50 {\~} # d # vert 50 {-} } ~ (U ~ stack { vert -50 
{\~} # d # vert 50 {-} } ) sup T]~ = ~U~E~[ stack { vert -50 {\~} # d # vert 50 

{-}}-( stack { vert -50 {\~} # d # vert 50 {-} } ) sup T ] - U sup T (10) 

By comparing equations (6) and (1 1 ), it can be seen that D is given by following 
20 equation: 

D ~ = ~ E ~ [ stack { vert -50 {V~} # d # vert 50 {-}}-( stack { vert -50 {\~} # d # 

vert 50 {-} } ) sup T ] 

25 

and that D, as defined above, is a diagonal matrix having 
the following format: 



30 
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D ~ = ~ SCALESYM 1200 [ MATRIX {{ stack { vert -50 {\~} # sigma }} & vert -50 
{0} & vert -50 {0) & vert -50 {...}& vert -50 {0} # vert -50 {0} & { stack { vert -50 
{\~} # sigma } } & vert -50 {0}& vert -50 {...}& vert -50 {0} # vert -50 {0} & vert - 
50 {0} & { stack { vert -50 {\~} # sigma }. } & vert -50 {...} & vert -50 {0} # STACK 
5 {.#.}& STACK {.#.}& STACK {.#.}& STACK {.#.}& STACK {.#.}# 
vert -50 {0} & vert -50 {0} & vert -50 {0} & vert -50 {...}&{ stack { vert -50 {\~} # 

sigma } } } SCALESYM 1200 ] 



Variances vert 50 {{ stack { vert -50 {\~} # sigma }} vert -50 stack {2#i}} are then 
determined from the diagonal matrix D. As can be seen from the above 
equations, the deviations in the vector { stack { vert -50 {Y-} # d # vert 50 {-} } } 
15 vert -30 sub i where i varies from 1 to N are uncorrelated and have the 
variances vert 50 {{ stack { vert -50 {\~} # sigma }} vert -50 stack {2#i}}. 

The final ^ value for determining a fault is obtained according 
to the following equation: 

20 d [ chi sup 2 ] ~ = ~ stack { vert -50 n # SCALESYM 200 SUM # {i=1}} {{{{ stack 
{ vert -50 {V} # d }} vert -50 stack {2#i}}}} over { vert 30{{ stack { vert -50 {\~} # 

sigma }} vert -50 stack {2#i}}} 



25 A normalization to = 1 as required in the definition of ^ 2 will be performed in 
equation (14). The value d[^ 2 ] is a single value which has reduced thermal 
and multipath noise, which represents information regarding a plurality of 
correlation measurements, and which may be compared to a threshold D in 
order to determine the existence of a fault. 

30 » Certain modifications of the present invention have been 

discussed above. Other modifications will occur to those practicing in the art of 
the present invention. For example, as described above, the ^ 2 distribution is 
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based on the assumption that all involved distributions are Gaussian. The 
distributions of d(< may deviate from this assumption and appropriate 

corrections to the formulas given here may be necessary. 

Moreover, the-present invention has been described above in 
5 connection with the detection of satellite signal faults such as those specified 
by the FAA, These faults result in signal distortions detectable by use of the 
present invention. The present invention as embodied by the following claims 
can also be used to detect other signal distortions such as those arising from 
multipath and satellite code cross correlation effects. 

10 Accordingly, the description of the present invention is to be 

construed as illustrative only and is for the purpose of teaching those skilled in 
the art the best mode of carrying out the invention. The details may be varied 
substantially without departing from the spirit of the invention, and the exclusive 
use of all modifications which are within the scope of the appended claims is 

15 reserved. 
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WHAT IS CLAIMED IS: 

The embodiments of the invention in which an exclusive property or right is 
claimed are defined as follows: 

5 1 . An apparatus for the detection of positioning system satellite 

signal distortions comprising: 

a correlator that determines a plurality of correlation 
measurements at points along a correlation curve, wherein each correlation 
measurement is based upon a correlation between a received satellite signal 
10 and a reference; and, 

a signal distortion detector that determines differences between 
the correlation measurements along the correlation curve and that detects a 
signal distortion from the differences. 

15 2. The apparatus of claim 1 wherein each of the correlation 

measurements represents a different time shift between the reference and the 
satellite signal. 

3. The apparatus of claim 2 wherein all of the different time 
20 shifts are late time shifts. 

4. The apparatus of claim 2 wherein all of the different time 
shifts are early time shifts. 

25 5. The apparatus of claim 2 wherein the different time shifts 

include a late time shift and an early time shift. 

6. The apparatus of claim 2 wherein the different time shifts 
also includes a zero time shift. 



30 
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7. The apparatus of claiin 1 wherein the signal distortion 
detector forms a deviation between the differences and expected values of the 
differences and compares the deviation to a threshold in order to detect 
existence of the signal distortion. 

5 

8. The apparatus of claim 1 wherein the signal distortion 
detector forms a deviation between each of the differences and a 
corresponding expected value of the difference and compares each of the 
deviations to a corresponding threshold in order to detect existence of the 

10 signal distortion. 

9. The apparatus of claim 1 wherein the signal distortion 
detector forms a deviation between each of the differences and a 
corresponding expected value of the difference, wherein the signal distortion 

15 detector determines a single deviation value based upon the deviations, and 
wherein the signal distortion detector compares the single deviation value to a 
threshold in order to detect existence of the signal distortion. 

10. The apparatus of claim 1 wherein the signal distortion 
20 detector forms a deviation between each of the differences and an expected 

value of the corresponding difference, wherein the signal distortion detector 
determines a covariance matrix based upon statistical properties of the 
deviations, and wherein the signal distortion detector uses the covariance 
matrix to perform a procedure on the deviations to create a single deviation 
25 indicative of the signal distortion. 

1 1 . The apparatus of claim 1 wherein the signal distortion 
detector detects the signal distortion by forming the following expressions: 

30 d sub {i,j} ~ - ~ Ed sub {i,j} 
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(Claim 1 1 continued) 

wherein each djj is a difference between a pair of correlation measurements i 
and j, and wherein Ed i? j is an expected value of the difference djj when there is 
no signal distortion. 

5 

12. The apparatus of claim 1 wherein the signal distortion 
detector detects the signal distortion in accordance with the following 
expressions: 

10 LINE d sub {ij} ~ - ~ Ed sub {ij} LINE ~ > ~ D sub {ij} 

wherein djj is a difference between a pair of correlation measurements i and j, 
wherein Edg is an expected value of the difference djj when there is no signal 
distortion, and wherein Dy is a threshold. 

15 

13. The apparatus of claim 12 wherein the signal distortion 
detector performs averaging in order to reduce effects of thermal and multipath 
noise. 

20 14. The apparatus of claim 12 wherein the signal distortion 

detector performs filtering in order to reduce effects of thermal and multipath 
noise. 



. 15. The apparatus of claim 1 wherein the signal distortion 
25 detector performs averaging in order to reduce effects of thermal and multipath 
noise. 

16. The apparatus of claim 1 wherein the signal distortion 
detector performs filtering in order to reduce effects of thermal and multipath 
30 noise. 
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17. A method of detecting signal distortions affecting a signal 
transmitted by a positioning system satellite comprising: 

correlating the transmitted signal with a first reference in order to 
determine a first correlation measurement at a first point along a correlation 
5 curve; 

correlating the transmitted signal with a second reference in 
order to determine a second correlation measurement at a second point along 
the correlation curve; 

correlating the transmitted signal with a third reference in order 
10 to determine a third correlation measurement at a third point along the 
correlation curve; 

determining a first difference from the first and second 
correlation measurements; 

determining a second difference from the second and third 
15 correlation measurements; 

directly comparing the first difference to a first threshold; 

directly comparing the second difference to a second threshold; 

and, 

detecting a signal distortion in the satellite based upon the 
20 comparisons of the first and second differences to the first and second 
thresholds. 

18. The method of claim 17 further comprising determining a 
third difference from the first and third correlation measurements and directly 

25 comparing the third difference to a third threshold, wherein the detection of a 
signal distortion comprises detecting a signal distortion in the satellite based 
upon the comparison of the first, second, and third differences to the first, 
second, and third thresholds. 
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19. The. method of claim 17 wherein the first, second, and third 
correlation measurements represent different time shifts between the reference 
and the transmitted signal. 

5 20. The method of claim 1 9 wherein all of the different time 

shifts are late time shifts. 

21 . The method of claim 1 9 wherein all of the different time 
shifts are early time shifts. 

10 

22. The method of claim 19 wherein the different time shifts 
include late and early time shifts. 



23. The method of claim 1 9 wherein the different time shifts 
15 also includes a zero time shift. 



24. The method of claim 17, 18, 19, 20, 21, 22, or 23 wherein 
the detection of the signal distortion comprises: 

forming a deviation between each of the first and second 
20 differences and a corresponding expected value for the difference; and, 

comparing each of the deviations to a corresponding threshold 
in order to detect existence of the signal distortion. 



25. The method of claim 17 wherein the detection of the signal 
25 distortion comprises: 

forming a deviation between each of the differences and a 
corresponding expected value of the difference; 

determining a single deviation value based upon the deviations; 

and, 

30 comparing the single deviation value to a threshold in order to 

detect existence of the signal distortion. 
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26. The method of claim 17 wherein the detection of the signal 
distortion comprises: 

forming a deviation between each of the differences and a 
5 corresponding expected value of the difference; 

determining a covariance matrix and mean values based upon 
statistical properties of the deviations; and, 

using the covariance matrix and mean values to perform a ^ 2 
procedure on the deviations to create a single deviation value indicative of the 
10 signal distortion. 

27. The method of claim 26 wherein the detection of the signal 
distortion comprises comparing the single deviation value to a threshold in 
order to detect existence of the signal distortion. 

15 

28. The method of claim 17, 18, 19, 20, 21 , 22, or 23 wherein 
the detection of the signal distortion comprises detecting the signal distortion in 
accordance with the following expression: 

20 d sub {i j} ~ ~ ~ Ed sub {i,j} 

wherein dy is the difference between correlation measurements i and j, and 
wherein Edy is the expected value of the difference djj when there is no signal 
distortion. 

25 

29. The method of claim 17, 18, 19, 20, 21, 22, or 23 wherein 
the detection of the signal distortion comprises detecting the signal distortion in 
accordance with the following expression: 



30 



LINE d sub {i,j} ~ - ~ Ed sub {ij} LINE - > - D sub {i,j} 



WO 02/39136 



PCT7US01/32163 



-21- 



(Claim 29 continued) 

wherein djj is the difference between correlation measurements i and j, wherein 
Edg is the expected value of the difference djj when there is no signal 
5 distortion, and wherein Djj is a threshold. 

30. The method of claim 17, 18, 19, 20, 21, 22, or 23 wherein 
the detection of the signal distortion comprises performing averaging in order to 
reduce effects of thermal and multipath noise. 

10 

31. The method of claim 17, 18, 19, 20, 21, 22, or 23 wherein 
the detection of the signal distortion comprises filtering in order to reduce 
effects of thermal and multipath noise. 

15 32. A method of detecting signal distortions affecting a signal 

transmitted by a positioning system satellite comprising: 

correlating the transmitted signal with references in order to 
determine a plurality of correlation measurements at corresponding points 
along a correlation curve; 
20 determining a single value from N values, wherein each value is 

formed based on a pair of correlation measurements, and wherein N > 2; 
comparing the single value to a threshold; and, 
detecting a signal distortion in the satellite based upon the 

comparison. 

25 

33. The method of claim 32 wherein each of the correlation 
measurements represents a different time shift between the references and the 
transmitted signal. 

30 34. The method of claim 33 wherein all of the time shifts are 

late time shifts. 
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35. The method of claim 33 wherein all of the time shifts are 
early time shifts. 

5 36. The method of claim 33 wherein at least one of the time 

shifts is a late time shift, and wherein at least one of the time shifts is an early 
time shift. 



37. The method of claim 33 wherein the different time shifts 
10 also includes a zero time shift. 



38. The method of claim 32 wherein the determination of the 
single value comprises: 

forming N differences between pairs of the correlation 
15 measurements; and, 

determining the single value from deviations between the N 
difference and corresponding expected values of the N differences. 

39. The method of claim 32 wherein the determination of the 
20 single value comprises: 

forming N differences between pairs of the correlation 
measurements; 

forming deviations between the N differences and expected 
values of the N differences; 
25 - determining a covariance matrix and mean values based upon 

statistical properties of the deviations; 

using the covariance matrix and mean values to decorrelate the 
deviations in order to form new deviations that are not correlated; and, 

performing a P> procedure on the decorrelated deviations to 
30 determine the single value. 
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40. The method of claim 39 wherein the use of the covariance 
matrix and mean values to decorrelate the deviations comprises using the 
covariance to form decorrelated and normalized deviations, and wherein the 

procedure is performed on the decorrelated and normalized deviations to 
5 determine the single value. 

41 . The method of claim 32 further comprising averaging in 
order to reduce effects of thermal and multipath noise. 

10 42. The method of claim 32 further comprising filtering in order 

to reduce effects of thermal and multipath noise. 

43. The method of claim 32 wherein the determination of the 
single value from N pairs of the correlation measurements comprises: 
15 defining a covariance matrix P in accordance with the following 

equation: 

P ~ = ~ E [ ( underline d ~ - ~ underline m ) ( underline d underline m ) sup 

T] 

20 

wherein the underlines indicate vectors, wherein E[A] is a statistical expectation 
of A, wherein the vector m is the mean value of the vector d, wherein the vector 
d is determined in accordance with the following equation: 

25 * {underline d} sup T ~ = ~ ( d sub 1 , d sub 2, d sub 3, d sub 4, . . d sub N) 

wherein N is the number of deviations, wherein the deviations d|< are formed 
from pairs of the correlation measurements lj and lj according to the following 
equation: 

30 

d sub K ~ = ~ I sub i ~ - ~ I sub j ~ - - Ed sub K 
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(Claim 43 continued) 

wherein Ed K is expected value of d K ; 

determining an upper triangular matrix U and a diagonal matrix 
5 D according to the following equation: 

p~ = ~UDUsupT 

defining stack { vert -50 {\~} # d # vert 50 {-} } in accordance with 
10 the following equation: 

stack { vert -50 {\~} # d # vert 50 {-}}- = - U sup {-1} ( underline d ~ - ~ 

underline m) 

wherein stack { vert -50 {\~} # d # vert 50 {-} } is a vector representing the 
1 5 decorrelated deviations generating the vector d; 

producing the following equation from the equations above: 

P~ = ~E~[U~ stack { vert -50 {\~} # d # vert 50 {-} } ~ (U - stack { vert -50 
{\~} # d # vert 50 {-} } ) sup T]~ = ~U~E~[ stack { vert -50 {\~} # d # vert 50 
20 B}~( stack { vert -50 {\~} # d # vert 50 {-} } ) sup T ] ~ U sup T 

determining the following equation from the equations above: 

D~ = ~E~[ stack { vert -50 {\~} # d # vert 50 {- } } ~ ( stack { vert -50 {\~} # d # 
25 vert 50 {-} } ) sup T ] 

wherein D has the following format: 



30 
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(Claim 43 continued) 

D - = - SCALESYM 1200 [ MATRIX {{ stack { vert -50 {\~} # sigma }} & vert -50 
{0} & vert -50 {0} & vert -50 {...}& vert -50 {0} # vert -50 {0} & { stack { vert -50 
{\~} # sigma } } & vert -50 {0}& vert -50 {...}& vert -50 {0} # vert -50 {0} & vert - 
5 50 {0} & { stack { vert -50 {\~} # sigma } } & vert -50 {...}& vert -50 {0} # STACK 
{.#.}& STACK {.#.}& STACK {.#.}& STACK {.#.}& STACK {.#.}# 
vert -50 {0} & vert -50 {0} & vert -50 {0} & vert -50 {...}&{ stack { vert -50 {\~} # 

sigma } } } SCALESYM 1200 ] 

10 

determining variances vert 50 {{ stack { vert -50 {\~} # sigma }} 
vert -50 stack {2#i}} from D above; 
15 determining a final ^ 2 value according to the following 

equation: 

d [ chi sup 2 ] ~ = ~ stack { n # SCALESYM 200 SUM # {i=1}} {{{{ stack { vert - 
20 50 {\~} # d }} vert -50 stack {2#i}}}} over { vert 30{{ stack { vert -50 {\~} # sigma }} 

vert -50 stack {2#i}}} 



and comparing d[^ 2 ] to a threshold D in order to determine 
25 existence of a signal distortion. 
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